Weak disorder strongly improves the selective enhancement of diffusion in a tilted 

periodic potential 
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The diffusion of an overdamped Brownian particle in a tilted periodic potential is known to 
exhibit a pronounced enhancement over the free thermal diffusion within a small interval of tilt- 
values. Here we show that weak disorder in the form of small, time-independent deviations from a 
strictly spatially periodic potential may further boost this diffusion peak by orders of magnitude. 
Our general theoretical predictions are in excellent agreement with experimental observations. 
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Diffusion plays a key role for mixing and homogeniza- 
tion but also for particle selection and separation tasks 
A particularly simple and common way to manipu- 
late the force-free diffusion of a Brownian particle is by 
means of a spatially periodic force field t2|| with a non- 
vanishing systematic component [1, 0, H, [al3i H H EIlj 
[T2I , i.e. the force derives from a tilted/biased periodic po- 
tential. Such dynamics arise in a large variety of different 
physical systems 0, H, H, Il0|. for example colloidal parti- 
cles in o ptic al potentials UQ, H, Q cold atoms in optical 
lattices 11 1, or globular DNA in microstructures p^ . 
While the force-free thermal diffusion of an overdamped 
Brownian particle is always reduced when switching on 
an unbiased periodic potential 0, the diffusion coeffi- 
cient as a function of an additional bias exhibits a pro- 
nounced peak 0, H, Q in a small vicinity of the so-called 
critical tilt, i.e. the threshold bias at which deterministic 
running solutions set in. This theoretical prediction has 
recently been confirmed by several experimental works 
01 HQ- However, for many experimental realizations of 
the above mentioned large variety of^stems involving a 
tilted periodic potential 0, i, i, i, 0, B i, UH, III , 
small, time-independent deviations of the potential from 
strict spatial periodicity are practically unavoidable. The 
objective of our present Letter is a detailed theoretical 
understanding of such weak disorder effects. An immedi- 
ate first guess is that they will somehow "wash out" the 
diffusion peak around the critical tilt. For instance, one 
might argue that enhanced diffusion requires a tilt close 
to criticality [9] and this fine-tuning will unavoidably be 
spoiled by the random variations superimposed to the 
original periodic potential. Here, we show that exactly 
the opposite is the case: Tiny deviations from spatial pe- 
riodicity result in an even more pronounced peak of the 
diffusion coefficient. Hence, the often unavoidable weak 
disorder is not an experimental nuisance but rather a new 
tool for sorting particles by way of a very strong and se- 
lective diffusion enhancement for certain species within 
a mixture. We remark that diffusion in the presence of 
temporal rather than spatial disorder represents a related 
but still different case. It also may result in accelerated 
diffusion, but, in contrast to our present case, already 



without a bias [Tj|. 

Our startingpoint is the usual overdamped Brownian 
motion in ID % i i, H, i, 0, i, S [3 : 



r]x(t) = -U'[x{t)) + y%kf^{t) , 



(1) 



where 77 is the viscous friction coefficient, and thermal 
fluctuations are modeled by unbiased, J-correlated Gaus- 
sian noise ^(t) with thermal energy kT. The potential 
U{x) consist of a tilted periodic part V{x) and "random" 
deviations W{x) (quenched disorder), 

U{x) = V{x) + W{x) , V{x) = Vo{x) ~ xF , 
Vo{x + L) = V^{x) , (2) 

where is a tilting force (static bias) and L the spatial 
period. Without loss of generality we focus on potentials 
"tilted to the right", i.e. F > 0. The quantities of main 
interest arc drift (average velocity) and diffusion. 
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where (•)^ indicates an average over the noise £,{t) in ([T]). 
While many of the following considerations can be gen- 
eralized to other types of disorder W{x), we focus on 
the analytically most convenient case of unbiased, homo- 
geneous Gaussian disorder. In other words, considering 
X as "time", W{x) is a stationary, Gaussian stochastic 
process with mean value {W{x)) — and correlation 



c{x-y) := {Wix)W{y)) 



(4) 



For simplicity only, we further assume that the correla- 
tion c{x) is monotonically decreasing for a; > from 



{Wixf) = c(0) 



(5) 



to c{x 00) = 0. A simple example is the critically 
damped harmonic oscillator 

W"{x) = -2A W'{x) - W{x) + 2ctA1/2 (g) 

driven by (5-correlated Gaussian noise "f{x), yielding a 
Gaussian W{x) with c{x) — (t^(1 + |a;|/A) e~l^l/'*', see 
Fig. 1. 
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FIG. 1: Diffusion D from ([Sjl versus tilt F for the dynamics 
©, © with Vo{x) = sin(a;) (i.e. L = 27r, = 1), 77 = 1, 
fcT = Dq = 0.001 (dimensionless units). The disorder W{x) is 
a stationary Gaussian process satisfying (O with cr = 0.03 and 
A = L. The insets illustrate (/(a;, F = 0) = Vo{x) + W{x) (ran- 
domness hardly visible) and W{x) (randomness alone) from 
(121). Filled dots (connected by short dashes): Numerically ex- 
act results. Dashed line: Analytical approximation ((Tj, (|15[). 
(fT7)) with (1231), (1211), and a = b = 0. Solid line: Exact analyti- 
cal result ((7ll-(l9l) in the absence of disorder (ct = 0). Lower and 
upper stars: analytical approximations (|22|l (no disorder) and 
(|25[) (with disorder), respectively. For large and small F (not 
shown) all curves monotonically approach D{F 00) — Do 
and D{F = 0) <C -Do. Main conclusion: the purely periodic 
potential leads to a maximal diffusion enhancement by about 
a factor of 70 (solid line) while a tiny amount of disorder 
(insets) further boosts the peak by a factor of 10 (dots). 



Without the disorder W{x), the following rigorous re- 
sults are known (see [0| and references therein) 

v = Do[l^ e-^^/'^^lM , D = Da B/A' , (7) 

where Dq :— kT/i] is the force-free diffusion coefficient 
according to Einstein and 
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dy dp dq 
/o L Jo Jo Jo 

V{x) - V{x ~y)- V{x -p) + V{x + q) , (9) 



see also 0, H, [13, for related findings. In particular, 
for F = one recovers the result 0] 



D{F = 0) 



Do 

C+C- 
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(10) 



Next we include the disorder W{x) 'm In a first 
step we assume that W{x + NL) = W{x) for some inte- 
ger N. Hence U{x) in ([21) is a tilted periodic potential 



with period NL. Accordingly, ((7l)- ((T0l) remain valid after 
replacing L by NL and V{x) by U{x). Exploiting 
one can rewrite (jSj after some manipulations as 

/''dye[^(-)-^(--«)l/'=^i(a:,y) (11) 



Jo L Jo 

A{x,y) ^(e[-«i^'+M^(-)->^(-f-"i)l/'=^^)(l2) 

n=0 

where {f{x)) ^~^E^L~oV(a; + i^L)- Next, we let 
N ^ 00 and adopt the usual tacit assumption [3, 
that this limit commutes with the limit t — > 00 in 
Exploiting that the Gaussian process W{x) is ergodic 



and satisfies (e ) = e 



{W^)/2 



yields 



A{X, y)=Y^ ^[-nLF+ciy+nL)]/kT (^3) 



Tl = 

,2 



c(a;) [ct^ - c(a;)]/fcr 



(14) 



From the monotonicity of c{x) (see above ([S])) one can 
infer upper and lower bounds for p3p and thus for (|11[) . 
Finally, one recovers the same formula for w as in ([71) but 
now with the relevant A given for any F > by 
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[V{x)^Vix^v)+c{y)]/kT 



(15) 
(16) 



with some (unknown) ■& £ [0, 1]. Note that a and c(y) are 
non-negative and that both quantities vanish (for all y) 
if and only if there is no disorder. In the latter case, (fT^l) 
reproduces ([51), otherwise the disorder always reduces the 
velocity v. This conclusion in fact applies to much more 
general types of disorder W{x), as can be inferred by 
applying Jensen's inequality to (fT^ . Finally, it follows 
from (fT5l) that v F/rj for F oo, independently of 
V{x) and W{x). For related findings see also p^ . 

An analogous calculation yields the same formula for 
the diffusion as in ^ with A given by (|15p and B by 

B=il + b) ^'^ ^dy ["^dp /'^dgels+'^l/'^-^ (17) 
Jo L Jq Jq Jq 

for any F > 0, where g is defined in ^ and 



h = c{y) + c{p) - c{q) - c{y-p) + c{y + q) + c{p+q) (18) 
6 - [1 - (1 - e-^^/'=^)3] [e-(-/fc^)' - 1] (19) 



with some (unknown) k S [—2,5]. One readily sees that 
1+6 > and that ^ is recovered in the absence of W{x). 
However, whether the disorder enhances or reduces the 
diffusion is not immediately obvious, with the following 
exceptions: From (fT^l) . ([TTI one finds that 

D{F ^oo) = Da (20) 

and from (fTUI - analogous to the derivation of (fT^l) - that 

D{F = 0) = (i?o/C+C_) e-('^/'=^)' . (21) 
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With Cauchy- Schwartz's inequality it follows flB*! that 
D{F = 0) < Do and with dTU]), (EB) that the disorder 
always reduces the dijfusion for F — 0. 

As exemplified in Fig. f and discussed in detail in 
0, without disorder, the diffusion D, considered as 
a function of the tilt F, develops a pronounced peak 
near Fc :— niaxj; V^'(x). The critically tilted potential 
V{x) — Vo{x) — xFc thus exhibits a strictly negative 
slope {V < 0) apart from plateaux {V = V" = Vq" ^ 0, 
V" — Vq" < 0) at Xc + nL for a generically unique 
Xc £ [0, L) and arbitrary integers n. In other words, in 
the noiseless dynamics ([l}, Fc marks the transition from 
locked to running deterministic solutions. For finite but 
weak noise, the peak of D{F) about Fc satisfies [§| 

D{Fc,W = 0) ~ 0.021 Do |Vo"'(xc)/fcTp/3 (22) 

With our findings D{F = 0) < Do and D{F ^ oo) ^ Do 
we thus can conclude that the peak height of D{F)/Do 
scales like T"^/'^, and similarly for its width 

Next we consider the influence of the disorder W{x) 
in the above most interesting regime of small F — Fc 
and kT. As intuitively expected and confirmed by closer 
inspection, in this regime the integrals in (|15p and (jl7[) 
are dominated by small values oi x — Xc y, p, q, thus 
admitting the following approximate expansions 

V{xc + 6) ~ V{xc)-{F-Fc)S + V^"{xc)5y6 (23) 
d{x) ~ c"(0)xV2 . (24) 

Further, we henceforth neglect a and b in and (fTT]) 
[l7| . The main virtues of these approximations are: (i) 
They still reproduce the correct limiting behavior for 
F oo and also qualitatively capture the small F be- 
havior, namely extremely small values of drift and dif- 
fusion. In other words, they are expected to reasonably 
work for all F > 0, as confirmed by Fig. I. (ii) The 
main effects of Vb(x) and W{x) are already captured by 
Fc, V^"{xc), and {W'{xY) = -c"(0). (iii) For F = Fc 
and sufficiently small kT, the integrands in and (fT7)) 
exhibit very pronounced maxima and thus can be evalu- 
ated by means of saddle point approximations, yielding 
the "universal scaling law" 

D{Fc) - D{Fc,W ^0) (1 + 1.9 Qe^^^'^'/^^ (25) 
Q := (I^'(a;)^)/|Fo"'(x,)p/3(fcr)4/3 . (26) 

The general analytical results \21\) . \25\) repre- 

sent the main findings of our present paper: Already a 
small amount of disorder typically leads to a much more 
pronounced peak of D{F) than without disorder. An 
illustration is provided by Fig. I. 

As an application of our general theory, we finally ad- 
dress the experiment from Ref. [1]: A colloidal sphere 
with diameter 1.48 /xm and -Do = kTroom/v — 0-19 /xm^/s 
moves along a ring of light. The particle feels = 80 
potential minima with period L ~ 0.33 /xm due to spa- 
tial variations of the light intensity and a torque due to 
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FIG. 2: Symbols: Experimentally measured diffusion D and 
velocity n (upper inset), adopted from Fig. 3 of Ref. 0]. Solid 
lines: Theoretical fit. Lower inset: Approximate intensity 
profile in arbitrary units for a — 1.0 and x £ [0, 20L]. 

orbital angular momentum transfer b y th e photons from 
two superimposed optical vortices [l,lll|, whose relative 
strengths are controlled by an experimental parameter 
a £ [0,1]. The resulting equation of motion takes the 
form ll]) ,18] with a potential (O corresponding to the 
total circumferential force from Eq. (4) in 4], namely 

- U'{x) = FoMa) + «-(«) cos(27ra;/L)] - W'{x) (27) 

with $(a) := (1 - a)/{l + a) and ^'(a) := 2[a{e'^^^ + 
C^)]^^V(1 + Further, Fq, e, and C 

are fit parame- 
ters, accounting for the laser intensity and the particle's 
shape, size, and composition (light scattering and ab- 
sorption properties). The torque being proportional to 
the light intensity gives rise to the first term on the right 
hand side of (P71) (fit parameter Fo) and also to part of the 
second term 18] (fit parameter e). Additionally, the sec- 
ond term accounts for the polarizable particle's coupling 
to the gradient of the light intensity (fit parameter (). 
Finally, W{x) in accounts for "random" imperfec- 
tions of the experimental optics [1] . Its variance is fixed 
by observingthat —W'{x)/Fo corresponds to the func- 
tion 77(6*) inQ] and that {\v{0)\'^) « 0.01 according to [1]. 
Further statistical properties of W{x) cannot be quanti- 
tatively related to those of the bare intensity reported in 
since both the intensity and the intensity gradient con- 
tribute - after suitably averaging over the particle volume 
- to W{x). For this reason, we model W{x) as Gaussian 
process © with periodicity W{x + NL) — W{x) and 
{W'{x)^) = (tVA^ = 0.01 Ftf. Regarding the correla- 
tion length A, we found that its exact quantitative value 
hardly matters and we have chosen A = 2L, in accordance 
with both Fig. Id from Q and the given particle size. 
Fig. 2 depicts our fit to the experimental results with 
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parameter values Fq = 1.37 pN, e = 0.38, and C = 0.25. 
[19| . In view of the periodicity W{x + NL) = W{x) we 
used the exact analytics (O with A from (fTTj) and simi- 
larly for B. The remaining dependence of the results on 
the realization of W{x) is still notable for N — 80. In 
the absence of more than one experimental realization, 
we have selected also in the theory a well fitting, but still 
representative single realization W{x). The minor differ- 
ences between theory and experiment can be naturally 
attributed to the fact that this W{x) is still not exactly 
the one realized in the experiment and to oversimplifica- 
tions of the theoretical model ([T]) per se. Once W(x) and 
all fit parameters are fixed in ([27)) . it is possible to ap- 
proximately estimate the underlying bare intensity I{x). 
The resulting I{x) in Fig. 2 indeed agrees quite well with 
Fig. Id from d. AU in all, our theory thus agrees in ev- 
ery respect very well with the experiment from Q. An 
analogous comparison with the experiments from [^, Q 
is prohibited by their small TV- values. 

In conclusion, our main finding consists in the general 
analytical results (PO)) . (PT|) . (^5)) . implying that even a 
tiny amount of disorder superimposed to the (dominat- 



ing) periodic potential may further boost the previously 
known sharp diffusion peak near the critical tilt by orders 
of magnitude, see Fig. 1. A further main point of our 
analytical findings is the universality of this very selec- 
tive and very strong diffusion enhancement close to the 
critical tilt. Considering that different species of parti- 
cles typically couple differently to the periodic and ran- 
dom potential and/or to the bias force we expect dif- 
ferent values of the critical tilt for each species. This 
opens the possibility of sorting particles by way of selec- 
tively enhancing the diffusion for certain species within 
a mixture. The experimentally often unavoidable weak 
disorder quite unexpectedly improves rather than deteri- 
orates the effectivity of the selection mechanism. Exper- 
iments along these lines for a mixture of different DNA- 
fragments in a periodically structured microfluidic device 
analogous to those from [13] are presently in preparation 
in the labs of D. Anselmetti at the University of Bielefeld. 
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